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ABSTRACT 


The BEocK version of the Splze=Step Fast Fourier 
Peanstocm (SSFFT) and the Jeager version of the Implicit 
Finite Difference (IFD) acoustic parabolic eguation models 
ace compared with a Lloyd mirror interference pattern in the 
mange domains, The SSFFT displays the inability to place the 
transmission loss nulis at the correct ranges. & 2S also 
unable to utilize bottom loss information correctly. The 
IFD produced nulls at the correct canges; however, IE 
inserted an unacceptable amount of noise except when small 
(1 m) vertical grid steps were used and the pressure release 
bottom was placed at extended depths. In shallow water 


cases, the IFD is able to properly represent the pressure 


information. Each model is explored in the wavenunmter 
domain cy use of a "Wavenumber Techniyue" (NT) model with 
emphasis on source depth determination. The source depth 


May be deternined by measuring the distance between the 
equally spaced nulls in the wavenumber representation. 
Neither acoustic model was able to provide accurate source 
depth information when the null spacings were compared toa 
knowh source-depth determination curve. Since the null 
spacings were not uniformly spaced, this was tc be expected. 
Some specific proplem areas in the modeis were identified by 
the use of the WT. 
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I. INTRODUCTION 
Conventional methods used to represent underwater sound 
transmission have deficted transmission loss (d3) as a func- 
eton of range in either tabular or graphical forn. Since 
the interest was in determining quantities such as range of 
bottom kounce, convergence zones, and “probabidaty of 
detection, these methods proved very satisfactory. However 
when it became necessary tc compare various underwater 
Beoustic propagation loss models with each other or actual 
in-situ experiments, they provided the analyst with less 
than adeguate insight into the causes for the observed 
errors in the various models outputs. Furthermore, interest 
has been generated in analyzing the received signal for 
information from the sound source which cannot easily be 
madauced from transmission loss as a function of range. 

One method, the Wavenumber Technigue (WT), described by 
Pens DiNapoli of NUSC [Ref. 1] and applied by Richard Lauer 
of NORDA [Eef. 2], provides more information for model 
comparisons and may have the capability of determining sound 
Source depth and range from the receiver in certain cases. 
DiNapoii used an analysis of acoustic propagation in the 
Wavenumber domain as an intermediate step in getting the 
EcanSmission loss curve in the Fast Field Program (FFP). 
Lauer studied the wavenumber domain information from the FFP 
and concluded that source localization might be vossible. 
Ee BS. Stamey [Ref. 3] investigated the potential for deter- 
Mining the source depth by using the Brock Split-Step Fast 
Fourier Transform (SSFFT) [Ref 4] parabolic eguation model 
with the "Wavenumber Technique" (WI) model. His preliminary 
investigation involved the use of isospeed and ASTREX sound 
fweed profiles with a fully absorbing botton, varied 


source/receiver combinations, and multiple freguencies. 
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The present investigation uses two acoustic parabolic 
eguation models for comparison with theoretical Lloyd mirror 
depictions, and the effects of observed inconsistencies in 
predicted transmission loss curves on the calcuiated results 
feom application of the WI model. An analysis of the sensi- 
tivity in determining source depth based on model output is 
attempted to ascertain possible use in naval operations. 

To facilitate a straight-forward analysis of the results 
and evaluation of the WT technigue for the cases studied, 
only a selected aumber of the various computer runs are 
included. Jnly those curves which best illustrate the 


effects revealed py this research are presented. 
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II, WAVENUMBER TECHNIQUE (WT) 


eee ae SS eee ee 


A. BACKGROUND (LLOYD MIRROR) 


The WI can be elucidated with the help of the classical 
Lloyd mirror effect [Ref. 5] which describes the interaction 
between direct path and surface reflected sound signals fron 
the same continuous wave (CW) source. Figure 2.1 illus- 
trates the geometry involved in the Lloyd mirror effect. 


ees 


Figure: 2.1 Lloyd Mirror Geometry. 
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The ccmplex pressure (P) can be expressed by 


1 -j z -j u 
psa| re ae eel ha 
: (eqn 2.1) 


when ©, and rz are related by 


R= ieee aees r= R= «(d+h)* (eqn 2.2) 


with the assumption r,~1r, the pressure amplitude (p) can be 


approximated by 


Ry = ain SUS ae (eqn 2.3) 


The interaction between the direct path and surface 
reflected waves produces constructive and destructive inter- 
ference which is manifested as peaks and nulls for the in 
phase and out of phase conditions respectively. Figure 2.2 
waS produced by using equation 2.3. The null locations can 


be obtained by the relationship. 


2) ee (ns | : 
pes ———— (egn 2.4) 
9) R sin R Yy 


When equation 2.3 is expressed in the horizontal wavenumber 
domain (k,), the details of the resultant functional depen- 


dency on horizontal wavenumber yield useful information. 
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B. PHYSICAL PROCESS 


The WI is a process by which the complex pressure wave 
(with range dependent amplitude and phase) is transformed 
mmeo the spectral density as a function of the horizontal 
wavenumber. The use of the WT in the operational environ- 
ment would require a quadrature demodulation of the source 
signal to attain the complex pressure. The quadrature 
demodulation process can best be defined by the illustration 


meetigure 2.3 [Ref. 6]. 


2cos (ut) 


Acos|(w -W)t+ x] 


LP 
filter 


Asin[(w ~G)t+ x| 


-2sin (axt) 
Figure 2.3 Quadrature Demodulation of a Signal. 
When acoustic models are used to simuiate the environment, 


the requirement to perform a quadrature demodulation is 


eliminated since the complex pressure is directly 
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accessible. Eo waGtidLae yy “aGoustwe parabolic, equation 
models provide the complex envelope function (U) which must 
be multiplied by the Hankel function in order to obtain the 
complex pressure. The acoustic wave iS written in the forn 


P=UHS(k rye ut (eqn 2.5) 


where the reference wavenumrer is defined by 


Ko= - (eqn 2.6) 
and 

BeUGs, Zz) (eqn 2.7) 
is the solution to the appropriate parabolic equation. The 


complex pressure is corrected for volume attenuation and 
Fourier transformed to attain the spectral density. A plot 
displaying the results graphically with the horizontal wave- 
Humber on the x-axis and the normalized spectral density on 
the yr-axis is constructed and the spacings between the nulls 
measured. Figure 2.4 {from Ref. 2] is an example of the NT 
for a source and receiver in the same type of water mass. 
Since the Lloyd mirror field in the wavenumber domain is 


given by [Ref. 2] 


E(k) = ae elie Ze) (eqn 2.8) 
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an alternative formulation [Ref. 2], which produces evenly 
spaced nulls, is obtained by converting the horizontal wave- 


number to beta 


BsAke = ke ails 


By using beta instead of tne horizontal wavenumber, the 
distance, delta beta, between any two nulls can be used to 


ascertain the source depth by 


maze (egn 2.10) 


as illustrated in figure 2.5. 

The application for implying tne method of images, in 
isospeed cases, to describe the erfects of a perfectly flat 
pressure release boundary can be justified by inspection. 
Therefore the parabolic equation can be used to produce the 
DPoyd's mirror effect and surface decoupling 1S not an 
issue, cf. Chapter 4 of Brekhovskikh [Ref. 7]. 

When the WT is implemented, care must be exercised to 
ensure that the complex pressure contains an adequate number 
Seepoints to describe the shortest periodicity in the pres- 
Sure ficid. Phas will preclude the possibility of aliasing 
in the wavenumber domain. The ceal and imaginary elements 
of the pressure signal must be modified so that the signal 
begins and ends with zero values. Piss WOdLeLeation 2's 
hecessary because the pressure signal must represent a 
repetitive oscillaticn for the transform. The complex pres- 
Pime array is zero-filled beyond the data in order to 
gemerate an array which has the length of a power of 2 for 
eae Fast Fourier Transform (FFT). The horizontal wavenumber 


increment is generated within the computer code by 
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III. ACOUSTIC MODEL ANALYSIS 


A. PARABOLIC EQUATICN 


The elliptical wave eguation can be approximated bya 
parabolic equation (PE) when it is assumed that the envelope 
function varies slowly with range. A detailed nathematical 
@escription of the PE from an acoustic point of view is 
presented by DeSanto [Ref. 8 and 9] and a simpler but less 
general description can be found in Coppens [Ref. 10}. The 
parabolic equation has the forn 


2 
= + 2)k (LY) + (k?-KE)U =O (ean 3-1 
Z 


where it is assumed that the pressure has the form of equa- 
elon 3.2, 


Pree 2) Sp) (eqn 3.2) 


S(t) represents the primary radial dependence of the field 
in terms of an outward propagating cylindrical wave of the 
form [Ref. 4] 


(act 
Se ao ules) (eqn 3.3) 


If it is assumed that the range of interest is many wave- 
lengths from the source; then the asymptotic form of the 
Hankel function (equation 3.4) can be used [Ref. 11]. 
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H61) 2 4/ eo e (kor -1/4) wikohns= 3 (egn 3.4) 


Two PE models, the Brock version of the Split-Sterp Fast 
Fourier Transform (SSFFT) and the Jeager version of the 
Implicit Finite Difference (IFD) were investigated for 
compatibility with the WT in determination of source depth. 


Be SSFFT 


1. Cansmission Loss Comparison with Lloyd 


LEEOL 


— == 


To analyze sound waves in the wavenumber domain, we 
needed a model to fprovide the pressure as a function of 
range. Initially the SSFFT model [Ref. 4] was selected . 
This model generates successive values for U as a function 
of range with the helr of the algoritha 

U(rear, 2) = e@Jarke(n?-1)/2 


(eqn 3.5) 


FRalech ojark? 
e JAPK2Ke FE T(U(r 7) 


This model was a natural choice since it was the model used 
in the preliminary study dealing with the determination of 
tae source depth by the WT [Ref. 3]. However, Zt Soon 
became apparent that this model's inherent weaknesses would 
require the consideration of another nodel if ocean botton 
interactions were to be studied. These weaknesses, and 
their impact upon the WT, will be discussed shortly. 

The SSFFT is a range dependent acoustic wave model 
which, for this analysis, will be operated in a range inde- 
pendent manner. in othe xs words, only one sound speed 
profile will be used, the water mass will be assumed homoge- 
neous, and the bottom flat. TO) facriitate Che study of the 
SSFFT model, the source code of the nodel was modified to 


generate the variables required for follow-on programs to 
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Seeamstorm the pressure information and then display the 
results graphically. LA copy Or tae, source code Listing for 
Gach program is provided in appendices A, B, and C. These 


programs were then linked by Job Control Language (JCL) so 
that the WT graphics were an automatically executed. 

A ieese rina t tne SSEPL was made tor 100 (nz, 500 mn 
source and receiver depths, 0.019 km range step, fully 
absorbing botton, 2000 m water depth, and an isospeed 
profile at 1500 n/sec. These values were chosen so that a 
comparison could be made with the theoretical Lioyd mirroc 
Beansmission loss (figure 2.2). Tre output was displayed as 
a@etransmission loss curve (figure 3.1). From equation 2.4, 
meemerccation o£ null number 1 (R1)’ should be 33.33 km and 
Hib: number 2 (R2) should be at 16.67 xn. The’ SshriLr placed 
fieadt approximately 38.3 km and R2 at 17.1 kno. Since the 
model kept the frequency, source depth, and receiver depth 
constant, then this equates to a reference sound speed of 
1305 and 1462 m/sec respectively. ft 2S “apparent fhatwa’s 
the range decreases the error is reduced which indicates the 
Bossibility of a problem in the range step section of the 
model which leads te the introduction of a systematic and 
increasing deviation of the calculated sound field from that 
Meedicted by the classical Lloyd's mirror interference 
pattern. Th2 “washing out" of the interference pattern for 
Bamges less tian about 2 km is a cesult of the particular 
appeoximation of a point source used to initiate the 
program, which is tc be expected. Ratner regular fluctua- 
tions develop beyond about 30 ka; from the geometry of the 
Sase, it is plausible that this arises from interference 
with bottom reflected signals. 
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Anotner test was made using the same inputs as above except 
that the source and receiver were nmaoved to a depth of 1500 2 
({tigure 3.2). Since the half beamwidth is 2) degrees and 
fme bottom is £ully absorbing, the shortest range at which a 
Peace reflection cculd occur is’ approximately 8.2 kn; 
beyond this range, surrace reflections should be spaced in 
accordance with equation 2.4. The nulls beyond 8.2 km do 
occur at the correct ranges; however, there are slow nodula- 
tions of the overall transmission loss signal. These siow 
modulations are another indication that the bottom is not 
mally aksorbing as had been specified as input to the 


progran. 


26 


J¢acta efor es # 
- wall to. ndaquh oe €S fave 
“vm 
tgs atpes (5 el rae” 
rs intew 9 089 2aad76ge oad ; 
-_ 
rewige ei Fees ib6> 


At 2608 sielgey ous Line 


UOTPOPISIUL woxgqog WOT TSpOW Laass po Soar, zee oanbty 


(WA) JONUY 
Ok AS 7h Sh ESecGme OPN ES 7p SC (| S2tz— az 


O3S/W O66°66hI O39dS GNNOS JONIYIIIY “990 000°02 


W/80, OX EZ2°Z 4900 NOILUNANSLLY “Wy 610°0 


STAI Sl Saal al Sez 


"YONI GINUY ‘HW 00° 0002 


ZH OO°OO! AONANOIY “W BI°OOST HLdIO YIAIIOIN ‘W OO*OOSI 


HLOIM WU48 41UH 


H1ddO 441M 
HLdI0 JIYNoOs 


(80) SSO7 NOISSIWSNYYL 


2) 


ee 


2. kT Analysis 


ine wi. Cutpue trom the SSFFT model 1s in Eng 
units, which is a probiem when comparing different nodels. 
However, for this investigation the problem was minimal. As 
expressed by equation 2.6, the reference wavenumber for 100 
hz with a reference sound speed of 4921.26 ft/sec (1500 
msec) is 0.1277 1/ft. After correction of the reference 
wavenumber for arrange increment of 62.3 it and the 2048 
points used to represent the wavenumber spectrun, ene 
Maximum beta value should be 0.1243 1/ft which agrees 
reasonably well with the maximum beta vaiues in figures 3.3, 
3.4, and 3.5. Contrary to theory the nuils are not equaily 
Seaced and the expected delta beta vaiue of 0.00192 1yit 
(Equation 2.10) is net represented by any of the null spac- 
ings. Stamey [Ref. 3], observed tne same non-uniforn 
Spacing of beta nulls. The reason why the source depth is 
not represented by any of the nuli spacings may be traced to 
the improper null placement in the raage domaiz. The Lack 
@£ egual spacing between the nuils must be attributed to 
something else, possibly the pressure wave was improperly 
represented as a continuous Signal or numerical errors exis 
W2Zthin the source ccde. Whatever the reason for this 


difference, it must Le resolved before the SSFFI can be used 
ie 


e 
epth. 


with the WI to determine the source dep 
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IFD 


the IFD nodel [Rref. Az 4 which was jnstalles and 
ssted inf September 4983 at NPS [ kef- 13] was the next nodel 

hosen as a sound source for the Wide This decision “ids made 
ecause the IFD was capable of properly handling the oceab 
orton qnteractions and was the only othet Acoust Te Vee nodel 
in cresidence at NPS- a copy of the source eoae Lasting cs 
provided rel appendix De 

The IFD also presented probiens which will be 
discussed further. The same initial test ran useé with the 
SSFFT was ysed with the IFD- Since the IFD requires more 
pottom jnfornation,s addit ional pesting we- nerformed in 
order to obtain 4 fui Ly apsorbing potton. Figure 3216 LS = 
trates the environment as seen DY the IFD model. 

The densities in the water maSS,y sediments and a6 ti- 
ficial attenuatiob layer were set to 4 constant of 1.0 and 
the attenuation in the sediment was adjusted until the elim 
qnation of bottom jnterfrerence was observed. Figure 324 2S 
an example of the IFD ¢transmissioa loss output for a 500 
source/receiver depth, #@4 water depth or 2000 He an uprer 
surface of the artificial attenuation layer at 3000 , and a 
lower pressure release surface at 1G00 Gs 22° bottom atten 


fieron used to produce figure 3.7 was 0.00 16 sa/wavelengt he 
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Figure 3.6 IFD Environment. 
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Mmmerocation of the nulls in figure 3.7 coincide with those 
predicted by equaticn 2.4 if the higa frequency noise is 
fapereda Out by either a running average process or low pass 
d@mgatal filter. A careful examination of the transmission 
Weegee curve (figure 3.8) reveals the real problem, which is 
the long wavelength noise near tae wavelength of interest. 
fme Long wavelength noise will cause erroneous spikes on the 
Meeplots after being Fourier transiormed. Continued anal— 
Mumemeecevealcd that -the noise could be eliminated if the 
thicknesses of the artificial attenuation layer and the 
bottom pressure release suriace are chosen so that they are 
at least twice the water depth. Phe, choice of a ‘small 
Memtical grid step will reduce the thickness reguired for 
Hee pressure release surface, but this version of the IFD 
Meegram 75 limited to 5000 vertical grid points due to soft- 
ware restrictions. FLgGUre 3.9 iS ,a transmission loss curve 
Obtained from a 500 m Source/receiver depth, a 2000 m water 
Bepth, a400) wm upper level artificial attenuation layer, 
and an 8000 m bottom pressure release surface. Comparison 
mmergures 3.1 and 3.9 with the classical Lioyd's mirror 


Peansmission loss (figure 2.2) shows that the shape of the 


curve beyond 30 km is greatly improved with the IFD. The 
Mumber of vertical grid steps was 5000. When the IFD is 
used in a sShallcw water situation and at 25 Z, the noise 


Mees not appear to be present. Figure 3.10 is a plot of the 
Beagsmission loss for a signal at 25 hz, with the sou 
Beceiver depth at 50 m, the bottom sloping up from 350 t 
M, and an artificial attenuation layer and pottom pressure 
Bebease surface located at 750 and 1000 om, respectively. 
Miesnoise which was so evident in tigure 3.7 is absent fron 
the shallow w«sater case. Producing tne noise free curves 
Tequires excessive computer time except for shallow water 


Cases. 
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To reduce the amount of the pressure information and 
Wee preserve the detaiis, the IFD tests were conducted 
between 10 and 40 km in range. The reference wavenumber for 
eae test input with the [FD was 0.4189 1/m which, after 
Correcting for a range increment of 7.5 m with 4096 pocints, 
produced a beta maximum of 0.4188 1/n. Pagures. S46 11, Ss. 12, 
and 3.13 indicate agreement with the beta maximum and a 
Metta beta of 0.00628 1/m can be found among the null spac- 
ings. When the vertical grid step is reduced in size, the 
Saeout Signal contains more information of finer detail. 
The finer detailed information will, if noise is present, 
Beagduce more clutter on the WT depiction (figure 3.13). in 
Numerous WI depictions a peak was observed at the nininun 
and maximum waveaumbers. Stamey believed that the right and 
left intensity maxima corresponded to beam elevation angles 
of 0 and 30 degrees, respectively. He further stated that 
the left maxima also may be reiated to algorithmic inaccura- 
eves [Ref. 3]. The "U-shaped phenomenon" observed during 
Seamey's investigaticn of the SSFFI, is present with the IFD 


and the nulls are not equally spaced. 
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fhe reason for the "U-shaped phenomenon" appears to be 
related to noise in the envelope function field which is 
further amplified when the Hankel function is inserted to 
obtain the conplex pressure. Using a low pass filter and 
increasing the vertical grid step size appear to reduce the 
'problen. tier giTe Seo, che mulls occur at 33.30 and 33.39 
| km men a Single null should exist at 33.33 km. This condi- 
tion or noise produces the most devastating erfects by 
| masking the desired density spectrun. A narrow band pass 
| filter would probably solve this problem in the short tern 
but eventually the IFD will require further study to elini- 
Mate the noise. When the shallow case is used in the WT 
(igure 3.14), the "U-shaped phenomenon" is absent and the 
mull locations are more distinct. However, in the shallow 
Water case, the bottom interaction produces high frequency 
anterference but the null Spacings occur near the expected 


Value for the source depth (0.0628 1/m). 
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Figure 3.14 


I¥. CONCLUSIONS 


iwo EalLabolic equation models, the SSFFT and IFD, were 
used to predict sound fields for comparison with the Llo 
Mirror interference fattern in the range domain and then in 
the wavenumber domain. In the range domain, th SSFFT 
mapEoperily handled the location of pressure nulls and 
displayed bottom interference. While these weaknesses in 
the SSFFI were absent in the IFD, the inadvertent and dele- 
terious insertion of noise was extensive. In tae wavenumber 
domain, the noise in the IFD was prohibitive in ail but the 
deep sediment and shallow water cases. The IFD produced 
results which were ccmparable with the SSFFI when the noise 
Was not a factor. The range step and the bottom attenuation 
aegorithms of the SSFFI should be reviewed in order to 
correct those problems described above. in the 25D, . the 
noise problem will require further investigation. ii tie 
problem in the iFD is not corrected, this model will be 
severely limited. With the guidance provided by Lauer 
[Ref. 2], the wavenumber domain can provide tne detailed 
information needed tO gusackiy correct inconsistencies 
observed between acoustic models. Therefore, the WT should 
be investigated further to ensure that any soitware weak- 
messes are eliminated and a thorough understazding of the 
process is assured. 

Experiments in the operational environment should be 
Planned to test the WT with actual sound sources. The ship- 
board experiment could be easily accomplished since the WT 
Can be iaplemented on micro-computers if the pressure infor- 
maewen 1S provided from an outside source. fhe WI, if it is 
going to be considered as a method to determine source 


depth, appears to have three major shortcomings: 
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1. A clean CW Signal with phase information is essen- 
tial, which will require quadrature demodulation and 
extensive filtering. 

2. Prediction models are required which will provide 
equally spaced beta nulis so that software can be 
developed to automate the process. 

3. High resolution processing equipment is necessary 
because the change in beta spacing for shallow sound 


sources is small (figure 2.5). 
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APPENDIX A: SSFFT SOURCE CODE (AT NPS) 
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